function model=r_model(dat,dataset,parameter,bbio,bbio_dir,year)
% ------------------------------------------------------
% -- Reading model output (following the observations)
%
% --- Datasets --- Kind --- Parameter ------------------
%    cruise        woce      T S NO3
%                  amt       T S CHL
%    mooring       bats      T S NO3 CHL PON DON PP
%    climato       woa05     T S NO3
%                  mld       MLD
%    satellite     seawifs   CHL
%                  eke       EKE
%                  sst       T
%                  pp        PP
%    floats
% ------------------------------------------------------
% G. Charria (02/2008)

switch dataset
   case 'cruise'
       [model]=r_nemo_cruise(dat,parameter,bbio,bbio_dir,year);
   case 'mooring'
       [model]=r_nemo_mooring(dat,parameter,bbio,bbio_dir,year);       
   case {'climato' 'satellite'}
       [model]=r_nemo_mx(dat,parameter,bbio,bbio_dir,year);       
end


